Chemically driven electron tunnelling pumps 
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The simplest mechanism for molecular electron pumps is discussed which is based on nonadia- 
batic electron tunnelling and nonequilibrium conformational fluctuations. Such fluctuations can be 
induced, e.g. by random binding of negatively charged ATP molecules to the electron-transferring 
molecular complex, their subsequent hydrolysis and the products dissociation. The pumping rate 
can be controlled by the ATP concentration in solution. Depending on the model parameters there 
may exist a critical ATP concentration for the pump to function. Alternatively, nonequilibrium 
fluctuations can be induced by externally applied stochastic electric fields. For realistically chosen 
parameters, the mechanism is shown to be robust and highly efficient. 

PACS numbers: 87.16.Uv,87.15.He,82.20.Gk,82.20.Uv,82.20.Xr,82.39.Jn,82.39.Rt 
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I. INTRODUCTION 



Electron transfer lies at heart of all bioenergetic pro- 
cesses. The energy of photoexcited electronic states, or 
one released in the oxidative breakdown of food molecules 
is used in a chain of electron-transfer reactions to create 
the transmembrane proton gradient storing ultimately 
the free energy ready to use in such "energetic" or- 

tnels of biological cells as mitochondria and chloroplasts 
1^ y, IJ, y{. This electrochemical proton gradient is 
used there by the ATP-syntase molecular complexes to 
synthesise the molecules of adenosintriphosphate (ATP) 
- a free-energy "currency" utilised in the most biochemi- 
cal cellular processes which occur far away from the ther- 
modynamic equilibrium under the living cell conditions. 
It is well known that the ATP-syntase can work also in 
reverse using the energy of ATP hydrolysis to restore the 
proton gradient |3]- The existence of the reverse elec- 
tron transfer, where the energy of ATP hydrolysis, or 
the free energy derived from electrochemical transmem- 
brane gradient of a sort of ions is used to energise the 
electroiis, is coming gradually in the focus of attention 
laSlalEllfl- The proton pumping molecular complexes 
driven normally by the energy released in the downhill 
electron transfer can work in reverse, pumping the elec- 
trons uphill on the time scale of seconds and minutes 
0, [l3- It might even be the case that such reverse 
electron transfer evolutionary emerged earlier in archae- 
bacteria existing at extremal environmental conditions 
(i.e., during the most earliest steps of the biological evo- 
lution). Moreover, the nitrogen fixation, which is realized 
by the nitrogenase protein complexes, utilises a p pare ntly 
the energy released in ATP hydrolysis [O, U2, HM Hjl • 
These natural molecular nanomachines produce ammo- 
nia routinely, at normal conditions, while the standard 
industrial technological process requires large pressures 
of about 150 atmospheres and temperatures in the range 
of 650-720 K. 

The electron transfer in nitrogenase provides one of 
the key steps in the overall reaction of ammonia syn- 
thesis. It is realized through a long-distance (about 14 
A) nonadiabatic electron tunnelling between two differ- 



ent metalloclusters situated in two different protein sub- 
units. This process is gated by nonequilibrium confor- 
mational transitions of the whole protein complex due 
to ATP binding and hydrolysis, with two ATP molecules 
hydrolysed per one electron transferred. How nitroge- 
nase works remains still a mystery, but a proper physical 
understanding gradually emerges [lj| . Such understand- 
ing is crucial not only for uncovering the working princi- 
ples of biological nanomachinery in general, but also for 
the molecular design of such and similar molecular ma- 
chines for a future nanobiotechnological use. It would al- 
low for an intelligent parameter optimisation, performing 
ultimately better than nature. Below I consider a very 
simplified, minimal theoretical model for such chemically 
driven electron tunnelling pumps which is based on the 
previous treatments in 14, 15, 16, 17, 18J . 



II. THEORETICAL MODEL 

A. Nonadiabatic electron transfer 

Let us start from nonadiabatic electron tunnelling cou- 
pled to molecular vibrational modes considered within 
a standard two-state, donor acceptor model (cf. one 
"frozen" conformation in Fig^J. This tunnelling can be 
described, e.g., within a spin-boson like model captured 
by the following Hamiltonian [H ES HI El HI IM El 

H = Ed\D){D\+Ea\A){A\ + Vtun{\D){A\ + \A){D\) 
+ li\D){D\-\A){A\)J2^,ib] + b,) 



-J2^u;,{byb, + l/2).il) 



Wherein, \D) and \A), are the localised electronic states 
with the energies Ed and Ea, correspondingly; Vtun 
denotes the effective electron tunnelling matrix ele- 
ment which incorporates the intervening medium influ- 
ence (a superexchange tunnelling mechanism is assumed 
El il El El EqI)- The coupling of electron tun- 
nelling to the molecular vibrational modes loj is char- 
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FIG. 1: Electron pumping scenarios based on dissipative elec- 
tron tunnelling and nonequilibrium conformational fluctua- 
tions. In the scheme I, the energy level of the donor state 
fluctuates in time. This is contrasted with the scheme II, 
where the acceptor level fluctuates in time. For dissipative 
tunnelling to occur from the donor site to the acceptor site of 
electron localisation, either the donor state should gain tem- 
porally in energy, or the acceptor cite should temporally lose 
in energy (conformation 2). A combination of two scenar- 
ios is also possible. The condition Kuni <C Vtun2 is required 
for the pumping mechanism to work robustly. It is necessary 
to block the reverse acceptor-to-donor electron transfer in the 
conformation 1. Temporal lifting of the electron energy is pos- 
sible, e.g. due to binding a negatively charged ATP molecule 
nearby the corresponding site of localisation. The tunnelling 
coupling can be exponentially reduced, e.g. due to increase of 
the tunnelling distance, or due to disruption of the tunnelling 
pathway. The latter one can be induced also by reorientation 
of a bridging molecular group. 



acterised by the coupling constants kj and the corre- 
sponding spectral density J{uj) — (27r/?i^) ^ k'jS{u! — 
ujj). This coupling modulates the energy difference 
between the donor and acceptor states. The corre- 
sponding fluctuations are described by (quantum) ran- 
dom force £,{t) = Y.j^^A^]^^''^'^ + ^je"'"^*) with ther- 
mally equilibrium autocorrelation function {^{i)S,{Q))T = 
n 

2 



The 



''— ^ J{uj)[coi\i[hijj/2kBT)coB{ijjt) — isa\ujt]dijj. 



medium's reorganisation energy X — h L dujJ{uj)/{2'KLj) 
serves is an integral characteristics of this coupling. An- 
other important medium's characteristic is the upper fre- 
quency ujc of low-frequency molecular vibrations, or sol- 



vent modes coupled to the electron transfer (ET), i.e. 
J{uj) = for a; > lOc- For T4un < A, y/kBTfiojc, the 
transfer kinetics occurs in the nonadiabatic tunnelling 
regime described by the quantum master equations of 
the Pauli type for the popula tions of donor and acceptor 

states miifliilillilMEi, 



Poit) = -kiPD + KpA, 
PA{t) = -kbPA + kf PD , 



(2) 



with the forward, donor-to-acceptor rate given by the 
quantum Golden Rule expression 



fcf = 



2V^ 



dTexp[-Q'(T)] cos [Q"(T)-er/n] , (3) 



where e — E^i — Ea is the difference of free energies 
(or the thermodynamic driving force —AG in chemical 
notations, e = —AG). Furthermore, the functions Q'{t) 
and Q"{t) in Q denote the real and imaginary parts of 
the doubly-integrated bath autocorrelation function with 
the reorganisation energy contribution added. 



Q(t) ^ dh m2)m)Tdh + iXt/h . (4) 

n Jo Jo 

The backward rate fch satisfies the Boltzmann relation 



fcb — fcf exp(— e/ZcsT) 



(5) 



for any J(w) |22l |. This yields equilibrium Boltzmann- 
Gibbs distribution p£)(oo)/pyi(oo) = exp(— e/fc^r) at 
temperature T, consistently with the condition of de- 
tailed balance, po{oo)kf — p/i(oo)fcb. It must be stressed, 
however, that the existence of the Boltzmann relations 
for rates like one in Eq. (0 at constant e and Vtun does 
not guarantee, generally, the detailed balance and the 
Boltzmann-Gibbs equilibrium distribution when these 
parameters explicitly fluctuate in time |25j . or when a 
stationary flux is present (given open, or cyclic bound- 
ary conditions even at constant rates). This is true even 
if the rates follow adiabatically to the instant values of 
the energy levels and the tunnelling coupling so that the 
condition (jSJ is valid at any instant of time. 

Independently of other details of electron-vibrational 
coupling, the quantum Golden Rule rate acquires in 
the high-temperature limit kgT ^ hiOc and the quasi- 
static approximation of electron energy fluctuations, 
{i{t)iiO))T « (I^(0))t « 2AfcRT the semiclassical 
Marcus-Levich-Dogonadze form [S^, IsJ 



fcf = 



2^ K? 



n ^AirXknT 



exph(e-A)V(4Afci3T)]. (6) 



Such a universality explains the widespread use of Eq. 
© in interpretation of experimental data. It involves 
(apart from temperature) three parameters only: the free 
energy difference e, the electron tunnelling coupling Vtun 
and the reorganisation energy A. When a coupling to 
high-frequency, essentially quantum modes is present, the 



exp ression © can be readily generalised accordingly [SJ, 
13^ . This presents one of the basic theories of electron 
transfer in molecular systems |23,|2J,|23. Of course, it is 
not truly universal, but it does provide a milestone and 
the simplest theoretical framework of practical relevance. 



B. Nonequilibrium conformational fluctuations 

Furthermore, let us assume that the electron- 
transferring protein complex (or its corresponding molec- 
ular subunit) can be in either of two conformations de- 
pending on binding a ligand, say ATP molecule (cf. 
Fig. ^. This assumption is quite in spirit of the 
Monod-Wwnan-Changeux model of allosteric enzymes 
|35l ISfl . ISTj . These two conformations possess very 
different I4un (distance between the donor and accep- 
tor sites is changed, or some bridging molecular group 
changes its orientation interrupting, or vice versa, estab- 
lishing thereby the electron-tunnelling pathway). These 
conformations can correspond also to very different en- 
ergy differences between the localised electron levels, e.g. 
the energy of the donor or acceptor state is changed 
(ATP and the hydrolysis products, ADP and the phos- 
phate group Pi, are all charged and the electrostatic ef- 
fects are of utmost importance here [ij). The attach- 
ment/detachment of ligand is a random process and both 
e{t) = Eoit) - EA{t) and Vtun(i) in Eq. QJ become 
stochastic functions of time. Alternatively, a strong two- 
state stochastic electric field can be externally applied 
to drive the electron transfer process. This is the start- 
ing point of the stochastically driven spin-boson model 
of Refs. [Ulie, 17, 18, 25j. It must be noted that mod- 
elling of the equilibrium conformational fluctuations in 
such a way should be considered with a great care [23 
(see also below), but the approach suits well to model 
non-equilibrium fluctuations like those considered |25l| . 
Within the approximations leading to Eqs. ©-Q, the 
quantum rates entering equations Q become stochastic 
functionals of Vtnnit) and e(t). Namely, the term er in 
Eq. (O is replaced by J^_^ e{t')dt' and instead of V^^^^ in 
the front of integral there appears Vtnn{t)Vtun{i—T) in the 
integrand [lj,[ig. However, if e(i) and Vtun(0 fluctuate 
slow on the time scale of Q(t), one can use an adiabatic 
driving approximation resulting in fluctuating rates fol- 
lowing to the instantaneous values of e(t) and Vtun(^)- 
This approximation is reasonable as a simple starting 
point for modelling and it can be justified in many cases. 
For these reasons, it is used below. The discussed adia- 
batic assumption means that after every conformational 
jump the vibrational relaxation to the new equilibrium 
of the vibrational degrees of freedom occurs very fast as 
compare with the mean duration of time spent in the 
corresponding conformation. Otherwise, the adiabatic 
driving approximation cannot be justified and the the- 
ory becomes essentially more intricate [ij, [1^, U^ [l3| • 

Furthermore, let us assume that the conformational 
fluctuations are Markovian and occur with the rates 



a and /? which do not depend on where the electron 
is localised, but are controlled by thermodynamically 
nonequilibrium concentrations of ATP, ADP, P; in the 
solution (i.e. ATP is continuously supphed |33|)- To 
be more concrete, let us assume that the conformational 
transition "1-^ 2" is caused by the ATP binding to 
the electron-transferring molecular complex (scheme I). 
Then, the transition rate a should obviously be propor- 
tional to the ATP concentration, [ATP], in the solution, 
i.e. a ex [ATP], since the binding frequency is propor- 
tional to [ATP] . On the contrary, the rate /3 of the con- 
formational transition "2— >1" caused by the ATP hydrol- 
ysis and the products dissociation should not depend on 
[ATP], but be rather determined by the activation bar- 
rier between two conformations and the energy released 
by breaking the phosphate bond. Such nonequilibrium 
fluctuations fuelled by this, or another source of chem- 
ical energy can drive electron transfer (FT) uphill. Al- 
ternatively, they can be induced by an externally applied 
stochastic electric field JS^, 113 • It can either be directly 
coupled to the electron transfer [l^ [13, 111], or modu- 
late the electron levels indirectly, via the electroconfor- 
mational coupling [391 ]. Then, within the discussed ap- 
proximations, the FT transfer kinetics is described by the 
kinetic equations Q with time-dependent rates undergo- 
ing two-state Markovian fluctuations. Formally, this is a 
typical problem of dynamical disorder j41i, [4^ 1431 144| . It 
can equivalently be described by the four-state Marko- 
vian kinetic scheme depicted in Fig. [21 ^45]. Similar 
schemes are standard by considering the problem of free 
energy transduction in biology [43, 113, IsJ • 
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FIG. 2: Equivalent kinetic scheme corresponding to quantum 
kinetic equations ([5J with rates undergoing two-state Marko- 
vian fluctuations 



Let us denote the population of the donor and acceptor 
states in the conformation i = 1, 2 as pDi and pAi, corre- 



spondingly. They obviously satisfy the master equations 



PDi = -(fcfi + a)pDi + PPD2 + KlPAl 

PD2 = apDl - iP + k{2)PD2 + kh2PA2 , 

PAl = kflPDl - (a + fcbl)pAl + f3pA2 , 

PA2 = fcf2P£)2 + apAl - (P + k\,2)PA2 ■ 



(7) 



These equations describe the process at thermodynam- 
ical equihbriuni if the overaU stationary flux is absent, 
i.e. clockwise and counterclockwise fluxes are mutually 
compensated. This requires that the product of forward 
rates along the cycle is equal to the product of backward 
rates, see e.g. in [S^, 111, |5J, |5J| . Otherwise, a nonequi- 
librium steady state (NESS) emerges with a persistent 
flux present. In such a case, one either requires a free 
energy supply to produce the corresponding stochastic 
cyclic motion (uphill motion on an effective free energy 
landscape), or this energy will be released (in the down- 
hill motion). Let us consider the situation where the 
averaged free energy bias e = eipf + £2^2* is negative, 
e < (pf* and P2* are the stationary probabilities of the 
corresponding conformations). Then the cycling in the 
counterclockwise direction in Fig. [21 is required to pump 
the electrons against the averaged free energy bias e < 
- the case of our interest here. If NESS corresponds to 
the clockwise total probability flux in Fig. |21 the roles of 
donor and acceptor states are interchanged and the con- 
formational fluctuations can be driven by the energy re- 
leased in the downhill electron transfer. These conforma- 
tional fluctuations can in turn be coupled to an ion flux 
to produce the uphill ion flow against the corresponding 
electrochemical gradient. This is the operating principle 
of the electron-driven proton pumps [3, Ull . The details 
are, of course, much more involved (a more complex, ex- 
tended kinetic scheme is required to describe these pro- 
cesses in a consistent manner) and still not completely 
understood. The operating principle is, however, rather 
clear due to (nonlinear) nonequilibrium thermodynam- 
ics considerations. In the present context, we combine 
them with a quantum treatment of the electron transfer 
kinetics. 

The free energy AGdrivc , which is required to drive one 
cycle on average in the counterclockwise direction in Fig. 
121 is determined by the well-known condition |5Q] 



fc^Tln 



AGdr 

fcb2fcfl 



fc^Tln 



knPkhia 



kh2f3hia 
= 62 — ei = Ed2 — Edi 



(8) 



For the scheme II, AGdrivc — Eai — Ea2- In the sec- 
ond line, we took Eq. jSJ into account. As discussed 
above, this energy can be delivered, e.g., due to the ATP 
hydrolysis (breaking the energy rich phosphate bond), 
or derived from any other free energy source (the pro- 
ton gradient, for example, via protonation/deprotonation 
of the molecular pump). Obviously, if only the tun- 
nelling matrix element is modulated by the conforma- 
tional transitions, then no pumping is possible and the 



described scheme agrees with the thermodynamic equi- 
librium, since no overall flux is present. We return back 
to the well-known problem of dynamical disorder in equi- 
librium systems and no more. To pump, one has to 
modulate the difference of electron energy levels e in 
time. With Vtun — const a pumping scenario conditioned 
on the existence of inverted Marcus regime of electron 
transfer, where the transfer rate decreases with the in- 
crease of the energy bias, is possible. It was described in 
Ref. [ial- In this respect, it is worth to notice that the 
very existence of the inverted FT regime presents a pro- 
foundly quantum-mechanical feature of nonadiabatic FT. 
Therefore, such an electron pump would be essentially 
quantum-mechanical. However, such a pumping scenario 
based on the sole modulation of e[t) would be rather in- 
efficient and too demanding for the system parameters 
in practice. Therefore, it is not likely to be used by na- 
ture. On the contrary, a properly concerted modulation 
of e{t) and Vtun(^) can allow to pump highly efficiently. 
In essence, for this one has to ensure Vtuni ^ Kun2 in 
Fig. ^ and a proper timing when FT kinetics is gated 
and locked to conformational fluctuations. Two possible 
pumping scenarios are depicted in Fig. ^ In the scheme 
I, the donor energy level is lifted upon binding negatively 
charged ATP molecule(s). Alternatively, one can modu- 
late the acceptor energy level in time, scheme II in Fig. 
n In the reality, a combination of both possibilities can 
take place. For example, in the case of nitrogenase the 
donor level is lifted by 300 meV and the acceptor level 
increases simultaneously by 100 meV, with the total in- 
crease of the driving energy bias by 200 meV |lj| ■ This 
compares well with the energy release from the hydrolysis 
of one ATP molecule which is about 0.3 — 0.5 eV under 
the living cell conditions. This basic pumping mechanism 
will be detailed and quantifled below. 



C. Solution of the model 

How to proceed further is standard and well-known 
|4a | . The solution of the master equations Q with the 
initial conditions Pj'(O) — 1, where f = Dl,D2,Al,A2 
yields the corresponding conditional probabilities of the 
state j, Pjj'{t). We are interested in several quanti- 
ties, such as (i) the asymptotic population of the donor 
state (pd(oo)); (ii) the time course of the donor state 
relaxation {po{t)), provided that the electron was ini- 
tially prepared in the donor state; (iii) the distribu- 
tion of the first arrival times, ipoiT), at the accep- 
tor state and the corresponding mean forward trans- 
fer time (t/) := J^ T'ipu{T)dT = J^ <I'£)(T)dT, where 
'^d{'t) = J il>D{T)dT is the corresponding survival prob- 
ability. Given the conditional probabilities Pjj'{t), the 
contracted probability of states Dl and D2 with the elec- 
tron being localised initially on the donor site is 



(pDlit)) = PDlDl{t)pf + PDlD2{t)pf 



(9) 



and 



{pD2{t)) = PD2Dl{t)pf + PD2D2(t)pf 



(10) 



correspondingly. In the above equations, it is tacitly as- 
sumed that the donor site has the same affinity to the 
transferring excess electron in the both protein confor- 
mations and these conformations are met with the sta- 
tionary probabihties pf = P/{a+(3) andp2* = a/{a+l3). 
The former assumption is trivially valid for the scheme 
II. However, for the scheme I it might be the case that 
the electron affinity to D2 state is much smaller than to 
Dl, i.e. the protein binds the transferring electron with a 
much higher probability in the first conformation (when 
no negatively charged ATP is bound nearby the donor 
state). In such a situation (which is not considered here 
for the sake of simplicity and analytical tractability of 
the results), one should put pf* -^ 1 and pf ^ in Eqs. 
© and p()|l . In any case, the averaged population of the 
donor state is {poit)) = {PDi{t)) + {pD2{t))- 

The formal solution can be found most conveniently 
using the Laplace-transform method. After some lengthy 
algebra we obtain (assuming equal electron binding affini- 
ties of the Dl and D2 states): 



(PDis)) 



lAis) 



sBis)' 



where 



A{s) = s^ + [v + pf {ki + ki,2) 



B{s) 



+pfik2 + kbi)]s 

+pfK2{ki + v)+pfk^i{k2 + v), 

s^ + {v + ki + k2)s + kik2 

+ (pffcl+pffc2)l/. 



(11) 



(12) 



(13) 



and i^ = a + /3, fci = fcfi -I- fcbi, fc2 — fcf2 + fcb2- In Eq. 
(|ll|l . {puis)) denotes the Laplace transform, (pd(s)) — 
/g exp(— st)(pD(t))di. The averaged asymptotic popu- 
lation of the donor level follows as (p£)(oo)) = A{0)/B{0), 



{pd{oo)) = 



P2*fcb2(fcl + v) +Pf fcbl(fc2 + v) 

kik2 + {pfki+p?^k2)v 



(14) 



The averaged relaxation of the donor state population is 
obtained by the inversion of Eq. Hll() to the time domain. 
It is bi-exponential and reads 

{pD{t)) = {poi^)) + [1 - (pD(oo))]i?(i), (15) 

where R{t) = Y^^^i 2 Cj exp(— F^t) is the relaxation func- 
tion with the rate constants 



Ti,: 



ki + k2 + ly 



± Viki + ' 



(3Y + AaP 



and the weighting coefficients 



1 



Cl,2 = 



1± 



ki+ k2 + a + P — 2co/ao 
^(fci +a-k2- PY+ AaP 



(16) 



(17) 



The remaining quantities oq and cq in Eq. (|17|l are 

Co = akf2{ki + a) + f3k{i{k2 + (3) + af3{kfi + k{2), 
Co ~ i''^{aki2 + (3k{i) + akikf2{2a + ki) 

+ Pk2kfii2l3 + k2)+a(3{ki+k2){kn + kf2). (18) 

The distribution of the first arrival times at the acceptor 
state ipoiT) can be immediately obtained from the sur- 
vival probability $£)(r) which in turn follows from the 
above relaxation function R{t) by setting fcbi7^b2 -^ 0, 
i.e. by assuming that the acceptor state is absorbing. 
This yield immediately 



$D(r) 



i=l,2 



Cjexp(-rir), 



(19) 



where the rate constants F^ and the coefficients Ci reduce 
to 



T,2='- 
2 



fcfi + kf2 + y 



± yj{kii + a-kf2- /3)2 + 4a(3 



and 



''^' = 2 



It 



a + (3+ (fcf2 - fcfi)(/3 - a)/{a + /?) 
y/{kn +a- kf2 - PY + 4a/3 



(20) 



,(21) 



respectively. These are the same expressions as, e.g., 
Eqs. (7)-(10) in Ref. ;48| obtained there using a different 
method. The corresponding mean forward transfer time 

is 



1^/; 



(a 4- I3f + akfi + pkf- 



(a-f /3)[afcf2+/3fcf 



knh 



'fl'Sf2 



(22) 



All the quantities, we are interested in, are thus formally 
determined. 



D. Conditions for pumping 

Let us suppose that the ET is characterised by Eq. 
© with the negative bias ei < and the tunnelling ma- 
trix element Vtuni in conformation 1 and the positive 
bias £2 > and the tunnelling matrix element Vtun2 in 
conformation 2. In addition, one assumes that the reor- 
ganisation energy A is the same in both conformations. 
Then for the ratio of the averaged donor and acceptor 
populations we obtain from Eq. H14|) : 



(Pg(oo)) 
{pa{oo)) 



exp 



£2 



l + g 
keTj l + C/i 



where 



pf 



{y + k2)F cos\\{e2/ 2k bT) 



pf z/cosh(e2/2fcBT) + k2Fcos\i{ei/2kBT) 
2 



(23) 



(24) 



F = ( 5^ ] exp[-(e? - el)/{A\kBT)] (25) 



Vt 



tun2 



TABLE I: Tunnelling coupling energies (in eV) and the cor- 
responding Marcus rates (in sec"^) 



Tunnelling coupling energies 



fcb 



Kunl = 5 ■ 10-% Kun2 = 10"* 
Kunl = 10-% Vtun2 = 10"* 
Ftunl = 5-10-% Kun2 = 5-10- 
Vtunl = 5-10-% Vtun2 = 2-10- 



0.54 29.5 1.18 ■ 10* 216 

0.02 1.18 1.18-10'' 216 

^ 0.54 29.5 2.95 ■ 10=^ 53.9 

5 0.54 29.5 471 8.63 



and ^ — exp[— (ei — e2)/(2fcBT)]. The pumping is most 
efficient when (p£)(oo))/(p^(oo)) <C 1. This requires 
£2 ^ ksT and (^ ^ 1. Moreover, to have the aver- 
aged energy gained by the transferring electrons maxi- 
mal, i.e. e = pfei + pfe2 ~ ei, one has to ensure that 
p|* ^ pf which contradicts, however, at the first look to 
the condition C ^ 1 in Eq. (|24|l . To resolve this contra- 
diction, one requires sufficiently small values of F ^ 1 
and this in turn demands 14uni ^ ^un2- How small is 
small depends on the system parameters. For example, 
for nitrogenase Ae = e2 — ei « 200 meV ,14]. Therefore, 
at the room temperatures, ksT « 25 meV, ^ « 54.6 is 
pretty large. Furthermore, let us assume for simplicity 
that |ei| = £2, so that C = iP/a){v + k2)F/{v + k2F) 
and F = (Vtuni/V^tun2)^- Moreover, we assume for a 
moment that a — 0.1/3, so that the conformation 1 is 
about ten times more probable than the conformation 2. 
Then, to satisfy ^^ <C 1 and to have an efficient pump- 
ing, Vtuni should be smaller than Vtun2 by, at least, two 
orders of magnitude. In such a case, fc2 ^ ^i, and for 
fci <^ :^ <^ ^2 , one can expect that the overall transfer will 
become gated by the conformation fluctuations. Namely, 
it follows from Eq. ^^ that for fcfi < a < /3 < fcf2, 
(r/) « 1/a [i^j, i.e. the pumping of electrons is locked 
to the conformational transitions caused by the binding 
of ATP molecules somewhere nearby the electron donor 
site (scheme I) .5^ , or vice versa by their hydrolysis and 
dissociation nearby the acceptor site (in the scheme II) . 
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FIG. 3: Dependence of the inverse of mean forward transfer 
time on the rate of conformational transitions a oc [ATP] 
(scheme I is assumed). 
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FIG. 4: Dependence of the ratio of donor and acceptor pop- 
ulations on the rate of conformational transitions a ex f ATPl . 



III. RESULTS AND DISCUSSION 

The outlined mechanism should be very robust. We 
illustrate it in Figs. El S for the following realistic test 
parameters which partially correspond to nitrogenase [l4| 
and partially are chosen just to demonstrate the essential 
effects: A = 1.2 eV, ei = -0.1 eV, £2 = 0.1 eV; the tun- 
nelling coupling energies are given in the TableQ]together 
with the corresponding Marcus rates. Furthermore, the 
value (3 = 1000 sec"^ is used in calculations and the rate 
a is varying as a control parameter assuming its propor- 
tionality to the ATP concentration [ATP]. 

As it is clearly seen in Figs. 13 0] the pumping 
effect is indeed present for Vtun2/Vtuni — 100 (con- 
tinuous line). Moreover, the mean forward time is 
locked to the rate a for a < 200 sec"^, being prac- 
tically its inverse. Moreover, the calculation of the 



survival probability ^^(t), e.g. for a — 100 sec ^, 
$£,(t) w 0.08 exp(- 12790 i) -f 0.92 exp(-92i); a = 200 
sec-\ ^D{t) ~ 0.14exp(-12798t) -I- 0.86exp(-184t), 
etc., shows that the transfer is almost single exponen- 
tial and (t/)-^ « a can be regarded as the pumping 
rate. With the further increase of a, the transfer kinetics 
becomes, however, ever more nonexponential with the ef- 
fective rate {Tf)~^ being larger than a, cf. Fig. O Here, 
the conformational transitions cease gradually to be the 
rate-limiting step and the tunnelling time 1/A:f2 becomes 
ever more important for the overall kinetics. However, 
this regime presents lesser interest in the present con- 
text since e = for a = f3 and the pumping effect then 
vanishes. 

Furthermore, let us to keep Vtun2 the same, but to 
increase I4uni such that the ratio I4un2/Vtuni becomes 
much smaller, Vtun2/Vtuni = 20 (fat dotted lines in Figs. 



130. This does not affect (t/) in Fig. 13 but changes 
dramatically the ratio of populations in Fig. 0] The 
pumping effect is present for sufficiently large rates a 
(cf. negative values of {p]:){(X))) / {pa{oo)) in Fig. 0J. The 
pumping efficiency drops, however, essentially. Moreover, 
the critical values of the rate ac and the associated ATP 
concentration [ATPc] emerge. The overall transfer oc- 
curs in the "donor— > acceptor" direction if only a > ac- 
On the other hand, one must keep a < (3. Otherwise, the 
transferred electrons will start to lose in energy on aver- 
age. Clearly, such a pump would not function perfectly. 
To realize a good electronic pump, the ratio Vtun2/Vtuni 
must be large. 

The inversion of the transfer direction depending on 
the rate of "1 -^ 2" conformational transition, cf. Fig. 
0] at several combinations of the tunnelling couplings, 
is rather intriguing. Namely, for sufficiently small a 
the roles of the donor and acceptor states are inter- 
changed. Here, the effective rate of the backward 
"acceptor— >donor" FT, defined as the inverse of the cor- 
responding mean first passage time (rt,), can be used to 
quantify the rate of transfer in this direction (assuming 
{pd{'X)))/{pa{'X))) ^ !)■ (Tb) can be obtained from Fq. 
(|22|l by setting there fcfi -^ fcbi,^f2 -^ fcb2- This quan- 
tity is depicted in Fig. [S| It is clearly seen in Fig. [S] 
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FIG. 5: Dependence of the inverse of mean backward transfer 
time on tlie rate of conformational transitions a. 



that the effective backward rate is not gated by the con- 
formational fluctuations for the used sets of parameters. 
Nevertheless, it can strongly depend on a, being almost 
linearly proportional to a at small a (see the continuous 
line in Fig. O. It might thus resemble a gating regime. 

One more interesting feature is that with the decrease 
of Vfau2 the effective transfer rate becomes smaller than 
a, with the tunnelling providing the rate-limiting step 
when a increases (see dashed line in Fig. |2l- 



A. Pumping efficiency 

The maximal pumping efficiency can be defined as the 
averaged energy gained by the transferred electron rel- 
ative to the energy required to drive one transfer cycle, 
i.e. ry = |e|/AGdrivc, or 



■q 



pf |ei| -p^*e2 



£2 



ei 



(26) 



For the above parameters (corresponding to the contin- 
uous lines in Figs. |31 El a-nd a ~ 100), the maximal 
pumping efficiency is rather high approaching 77 = 0.5 
56]. It can be even higher approaching one, if the affin- 
ity of the donor state Dl to electrons is much larger than 
the affinity of the state D2 (scheme I), i.e. the protein 
complex takes preferably electrons from the bulk in the 
conformation I (the formal solution of the model has to 
be modified in this case, but the qualitative features re- 
main). 



IV. CONCLUSIONS 

The considered generic model might seem somewhat 
oversimplified. It is indeed aimed primarily to highlight 
the basic working principles and their practical relevance. 
This model should be extended and generalised further in 
several directions, e.g. a correlated two-electron transfer 
should probably be considered in nitrogenase as an el- 
ementary step rather than single-electron transfer and 
a proper treatment of the ATP binding, hydrolysis and 
dissociation of the hydrolysis products would require to 
introduce more conformations than two. Moreover, the 
external uptake and release of electrons from and to the 
donor and acceptor sites, e.g., from mobile electron carri- 
ers should be incorporated in the complete model. Nev- 
ertheless, the considered elementary model does allow 
to manifest the main operating principles which are not 
much different from those well established and clearly 
understood, both phcnomenologically and in progressing 
details, for ionic pumps 0j|j)l2y|- Moreover, it allows 
one to clarify some important conditions for the efficient 
pumping such as a large ratio of the tunnelling couplings 
in the different conformations and a possible existence 
of the critical ATP concentrations. The profound physi- 
cal difference between the ionic and electronic pumps is, 
however, that the electron is essentially a quantum par- 
ticle and it tunnels over a large distance between met- 
alloclusters in nitrogenase (also in other electron trans- 
fer complexes, like cytochrome bci) using virtually pro- 
tein bridging states. This is why the details here are 
definitely very different from ionic pumps. They do 
matter and are important to arrive in a future at the 
detailed (quantum)-mechanistic, molecular-dynamic un- 
derstanding which still is lacking at present. Unlike to 
the many- years, extensive research on ionic pumps we 
undertake here really the first steps. The research do- 
main of electron transfer driven by a chemical energy 



source through noncquihbrium conformational fluctua- tions is just emerging. 
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